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A new physical mechanism of formation of runaway electron (RE) beams during plasma disrup¬ 
tions in tokamaks is proposed. The plasma disruption is caused by a strong stochastic magnetic 
field formed due to nonlinearly excited low-mode number magnetohydrodynamic (MHD) modes. 

It is conjectured that the runaway electron beam is formed in the central plasma region confined 
inside the intact magnetic surface located between q — 1 and the closest low-order rational magnetic 
surfaces [q = 5/4 or q = 4/3, ...]. It results in that runaway electron beam current has a helical 
nature with a predominant m/n = 1/1 component. The thermal quench and current quench times 
are estimated using the collisional models for electron diffusion and ambipolar particle transport in 
a stochastic magnetic field, respectively. Possible mechanisms for the decay of the runaway electron 
current owing to an outward drift electron orbits and resonance interaction of high-energy electrons 
with the m/n = 1/1 MHD mode are discussed. 


The runaway electrons (REs) generated during the dis¬ 
ruptions of tokamak plasmas may reach a several tens of 
MeV and may contribute to the significant part of post¬ 
disruption plasma current. The prevention of such RE 
beams is of a paramount importance in future tokamaks, 
especially in the ITER operation, since it may severely 
damage a device wall [l|-jf|. 

The mitigation of REs by massive gas injections (MGI) 
and externally applied resonant magnetic perturbations 
(RMPs) have been extensively discussed in literature 
(see, e.g., Refs. [fUU and references therein). However, 
no regular strategy to solve this problem has been de¬ 
veloped because up to now the physical mechanisms of 
the formation of REs during plasma disruptions are not 
well understood. In spite of the numerous dedicated ex¬ 
periments to study the problem of runaway current gen¬ 
eration during plasma disruptions in different tokamaks 
(see, e.g., |7Hla|) no clear dependence of RE formation 
on plasma parameters has been established. These nu¬ 
merous experiments show the complex nature of plasma 
disruption process especially the formation of RE beams. 

One of the important features of the formation of RE 
beams is the irregularity and variability of the beam pa¬ 
rameters from one discharge to another one. This is 
an indication of the sensitivity of RE beam formations 
on initial conditions which is the characteristic feature 
of nonlinear processes, particularly, the chaotic system. 
Therefore one expects that ab initio numerical simula¬ 
tions of the RE formation process may not be quite pro¬ 
ductive to explain it because of complexity of computer 
simulations of nonlinear processes [161 ] . The problems of 
numerical simulations of plasma disruptions is compre¬ 
hensively discussed in [la¬ 
in this work we propose a new physical mechanism 
of formation of RE beams during plasma disruptions in 
tokamaks. It is based on the analysis of numerous exper¬ 
imental results, mainly obtained in the TEXTOR toka¬ 


mak and the ideas of magnetic field stochasticity jl8[ . 
The mechanism explains many features of plasma dis¬ 
ruptions accompanied by RE generations. 

Main conjectures. It is believed that the plasma 
disruption starts with the excitation of magnetohydrody¬ 
namic (MHD) modes with low poloidal m and toroidal n 
numbers, (m/n = 1/1, 2/1,3/2, 5/2,...) that lead to a 
large-scale magnetic stochasticity (see, e.g., (Tol [22[ and 
references therein). The heat and particle transports in 
the strongly chaotic magnetic field causes the fast tem¬ 
perature drop and ceases the plasma current. This pro¬ 
cess depends on the structure of the stochastic magnetic 
field which depends on the spectra of magnetic pertur¬ 
bations and on the safety factor profile q(p), (p is the 
minor radius of the magnetic surface). At certain condi¬ 
tions the stochastic magnetic field may not extend up to 
the central plasma region due to the creation of the outer¬ 
most intact magnetic surface p c . The electrons confined 
by this magnetic surface are accelerated by the toroidal 
electric field induced by the current decay from the outer 
plasma region, which leads to the formation of the RE 
beam. The initial RE current Ip RE ^ is mainly determined 
by the pre-disruption plasma current distribution I p (p) 
confined by the outermost intact magnetic surface p c , i.e., 

4 RE) « i p (pc). 

The lifetime of the RE beam mainly depends on two 
effects: the outward drift of RE orbits induced by the 
toroidal electric field E v [H, 33] and the resonant inter¬ 
actions of REs with helical magnetic perturbations. The 
first one is responsible for the smooth decay of the RE 
current, while the second one cause the sudden RE losses. 
The outward drift velocity Vdr is determined by E v and 
the RE current, Vdr oc E v> /Ip RE ' > oc E v / p 2 [32), [33j. The 
most stable of the RE beams are expected to form when 
the corresponding drift velocity is lowest and the low- 
order rational surfaces within the RE beam are absent or 
one. 
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Consider, for example, the pre-disruption plasma with 
a monotonic safety factor profile q(p) with q( 0) < 1. 
Then the most stable RE beam can be formed when 
the outermost intact magnetic surface is located between 
magnetic surface q = 1 and the nearest low-order ra¬ 
tional surfaces q = 5/4 [or q = 4/3, ...]. It occurs at 
the sufficiently small amplitude of the m/n = 1/1 mode. 
There is only one rational magnetic surface q = 1 within 
the RE beam that is resonant to the large-scale mag¬ 
netic perturbations, particularly, to the RMPs. Such RE 
beams are relatively stable since low-energetic REs (up 
to 10-15 MeV) are not destabilized due to absence of a 
large scale stochasticity. The loss of REs mainly occurs 
due to the outward drift of RE orbits and the stochastic 
instability of high-energetic REs due to the interactions 
of high-mode harmonics of the m/n = 1/1 mode of mag¬ 
netic perturbations. 

In the case of plasma disruptions with q( 0) > 1 the 
intact magnetic surface p c would be smaller while the 
toroidal electric field E v would be larger than in the ones 
with g(0) < 1. Due to the large outward drift velocity 
Vdr such RE beams would cease faster. 

The two possible distinct generic structures of a 
stochastic magnetic field before the current quench with 
the RE-free discharge and with the RE discharge are 
shown in Figs. [T| (a) and (b) by the Poincare sections 
of magnetic field lines. It is assumed that the pertur¬ 
bation magnetic field contains several low-mode number 
m/n MHD modes with equal amplitudes B mn : (a) the 
amplitude Bn of the m/n = 1/1 mode is equal to oth¬ 
ers; (b) B n is four times smaller than the amplitudes 
of other modes. As seen from Fig. [T| (a) for the large 
amplitude of the (m/n = 1/1) mode the stochastic mag¬ 
netic field extends up to the central plasma region de¬ 
stroying the separatrix of the m = n = 1 island. For 
the low-amplitude of the ( m/n = 1/1) mode shown in 
Fig. [T] (b) the stochastic magnetic field does not reach 
the q = 1 magnetic surface. The last intact drift surface 
(red dots) is located between the resonant surfaces < 7=1 
and q = 5/4 (blue curves). 

The existence of an intact magnetic surface and its lo¬ 
cation depends on the radial profile of the safety factor 
and on the spectrum of magnetic perturbations. The 
latter sensitively depend on the plasma disruption condi¬ 
tions and vary unpredictably from one discharge to an¬ 
other during plasma disruptions. This makes RE forma¬ 
tion process unpredictable and may explain a shot-to- 
shot variability of the parameters of RE beams. 

Experimental evidences. This conjecture on the 
mechanism of RE beam formation agrees with the im¬ 
portant features of the experimental observations in the 
TEXTOR tokamak. In the experiments the plasma dis¬ 
ruptions were triggered by gas injections (see, e.g., J|- 
UlL 1231 ): the disruptions with REs were triggered by ar¬ 
gon (Ar) injection and the RE-free disruptions with Ne 
injection. The injection of these gases may finally give 
rise to different spectra of amplitudes of MHD modes. 
One can expect that the amplitude of the m/n =1/1 
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FIG. 1: (color online) Poincare sections of magnetic field 
lines in a pre-disruption plasma caused by several m/n MHD 
modes, (n = 1,2, 3; m = 1, ...8): (a) all mode amplitudes 
Bmn are equal; (b) the amplitude Bn of the m/n = 1/1 
mode is four times smaller than Bmn■ The safety factor at the 
magnetic axis q( 0) = 0.8 and at the plasma edge q(a ) = 4.7. 


MHD mode excited by the He/Ne injection is higher than 
in the case of Ar gas injection. 

The plasma current decay in the current quench (CQ) 
and the RE plateau regimes for all discharges is well ap¬ 
proximated by the linear function of time I p = I p o + bt , 
with the average CQ rate b = (dl p /dt ) as shown in Fig. [2] 
(a). The current decay rates \{dl p /dt)\ in the CQ stage 
and the RE plateau stage versus the initial RE current 
I p RE ' > for a number discharges are plotted Fig.[2](b). The 
plausible radial profiles of I p (p) and the corresponding 
safety factor q(p) are plotted for the two values of g(0) 
in Fig. H 

Existence of the finite interval of the initial RE cur¬ 
rents I p RE \ Since p c is located between the magnetic 
surfaces pi and p$ corresponding to q(pi) = 1 and 
q(ps) =4/3, the RE current I p ' should take values in 
the finite interval. This expectation is supported by the 
experimental values of the plasma current l p RE> as seen 
from Fig. [2] (a) and (b). These values of I p RE ' > also lie 
in the region between the resonance magnetic surfaces 
q(pi) = 1 and q{pz) = 4/3 [or g(p 2 ) = 3/2] as shown 
in Fig. [3] where the radial profile of the pre-disruption 
equilibrium plasma current I p (p) (curve 1) and the cor¬ 
responding safety factor profile q{p) (curve 2) are plotted. 

The average values of {\dl p /dt\) for almost all dis¬ 
charges are confined in the interval (2.2, 5.6) MA/s, i.e., 
in one order lower than the current decay rate in the CQ 
stage. The values of I P RE are in the range between 170 
kA and 260 kA (see Fig.[2](b)). These values of {\dl p /dt\) 

and I p RE ' > are close to the ones observed in the similar 
experiments in the DIII-D tokamak (see, e.g., [l5jl. 

As seen from Fig. [2] (a) there are untypical discharges 
with the highest and lowest values of I p RE ’ that corre¬ 
spond to p c at the borders of region pi < p < p 3 . For 
these discharges the CQ rates (\dl p /dt\) take highest or 
lowest values. The RE current decay rates of these dis- 
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FIG. 2: (color online) (a) Typical time evolution of the plasma 
current with RE current. The average current decay rates 
(dl p /dt) at the CQ and the RE plateau stages are determined 
by fitting with a linear function I p (t) = a + bt. Symbol 0 
corresponds to the plasma current I P RE ^ at the initial stage of 
the RE plateau, (b) The decay rates \{dl p /dt)\ versus I p RE \ 
Symbols □ correspond to the CQ rate (l.h.s. axis), and © — 
the RE plateau (r.h.s. axis). 


charges take the highest values. They have the shortest 
duration time of RE currents. One expects that the pres¬ 
ence of several low-order m/n = 4/3, m/n = 3/2, and 
m/n = 1/1 resonant magnetic surfaces within the RE 
beam for the discharge with the highest Ip RE ^ may lead 
to excitations of the corresponding MHD modes. The 
interactions of these modes may lead to the quick loss of 
REs due to the formation of a stochastic zone at the edge 
of the RE beam. 

Dependence on the level of magnetic perturbations. 
The existence of the intact magnetic surface p c between 
the q = 1 and q = 4/3 rational magnetic surfaces and 
its location depends on the level magnetic perturbation 
cmhd (more exactly on the spectrum B mn ). With in¬ 
crease of cmhd the radius p c shrinks and it can be bro¬ 
ken at the certain critical perturbation level e cr . It leads 
to the total destruction of confinement of electrons and 
ions. This is in agreement with experimental observa¬ 
tions of the existence of critical magnetic perturbations 
from which on runaway beams are not generated 0 - 

The shrinkage of p c with increasing the magnetic per¬ 
turbation cmhd leads to the decrease of the RE cur¬ 
rent Ip RE ' > since I p RE « I p (p c ). On the other hand if 
one assumes that the plasma current decay is caused by 
the radial transport of particles in the stochastic mag- 



FIG. 3: (color online) Radial profile of the plasma current 
I p (p) (solid curves 1 on l.h.s. axis) and the corresponding 
safety factor profile q(p) (dashed curves 2 on r.h.s. axis). 
The rectangular (red) dots correspond to the experimentally 
measured values of I P RE ^ for several TEXTOR discharges. 
The plasma parameters are I p = 350 kA, Bo = 2.4 T, Ro = 
1.75 m, and a = 0.46 m. The values of q(0) are 0.7 (solid black 
curves) and 0.8 (dashed magenta curves), respectively. The 
radii pi, p 2 , and po are the positions of the rational magnetic 
surfaces q(pi)=l, q(p2)=3/2, and q(po)=A/3, respectively. 


netic field, the CQ rate dl p /dt should be proportional 
to the square of the magnetic perturbation level cmhd-, 
\(dl p /dt}\ oc \cmhd\ ■ Therefore, one expects that to the 
higher values of | ( dl p /dt) | correspond the lower values of 
the RE current I p RE K This expectation is in agreement 
with the experimental measured values of these quanti¬ 
ties presented in Fig. [2] (b). 

The formation of the RE beam inside the intact mag¬ 
netic surface can be also confirmed by the spatial profiles 
of the synchrotron radiation of high-energy REs with 
energies exceeding 25 MeV. One observes that the radia¬ 
tion is localized within a finite radial extent in the central 
plasma region. 

Estimations of thermal quench (TQ) and CQ 
times. The strong radial transport along the stochastic 
magnetic held lines causes the losses of heat and plasma 
particles from the stochastic zone. The TQ can be ex¬ 
plained by the fact that the anomalously large heat trans¬ 
port in a stochastic magnetic held is mainly determined 
by the electron diffusion. The CQ is determined by the 
particle transport in a stochastic magnetic held and has 
an ambipolar nature. Using the collisional test particle 
transport model in a stochastic magnetic held m we es¬ 
timated the heat conductivity Xr(p) and the ambipolar 
diffusion coefficient D p of particles. 

For typical magnetic perturbations and pre-disruption 
plasma temperatures (0.5 4- 1.0 keV) the magnitude of 
Xr{p) has the order of several 10 2 m 2 /s. The charac¬ 
teristic heat diffusion time th = a 2 /2yy is of the order 
of 10 -4 s that agrees with the experimentally observed 
times. The quantitative analysis based on the numerical 
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solution of the heat diffusion equation also gives similar 
values for r#. 

The ambipolar particle transport in a stochastic mag¬ 
netic field is strongly collisional due to the low plasma 
temperature (from 5 eV to 50 eV) after the TQ. At 
these plasma temperatures the corresponding diffusion 
time t p = a 2 /D p of particles changes from 1 s to 0.3 
s. Since the diffusion coefficient D p oc B/ nn and there¬ 
fore t p oc B “ 2 , then t p can be reduced to one order 
smaller value for a three times larger perturbation than in 
Fig.m This timescale is still much longer than the exper¬ 
imental values. However, the collisional model does not 
takes into account the effect of the toroidal electric field. 
One expects that the acceleration of electrons and ions 
by the toroidal electric field increases the radial trans¬ 
port of particles. To include this effect in the collisional 
model one can assume that the effective temperature of 
the plasma is higher than the measured one. The parti¬ 
cle diffusion time t p at the effective temperature 2 keV is 
about 8 x 10~ 3 s. This timescale gives the average current 
decay rate d,I p /dt « I p /t p = 0.35/(8.0 x 10 -3 ) « 44.0 
MA/s which is order of the experimental measured one 
given in Figj2](b). 

In general the transport of heat and particles in the 
presence of RMPs is a three-dimensional problem. Par¬ 
ticularly, a stochastic magnetic field with the topological 
structures like the ones in Figs. [T| leads to poloidally and 
toroidally localized heat and particle deposition patterns 
on the wall (see, e.g., 0 ) similar to those in ergodic 
divertor tokamaks (see, e.g., 0 )- 

RE beam evolution. From the described scenario 
of plasma disruption it follows that a typical runaway 
beam current is localized inside the area enclosed by the 
last intact magnetic surface. In general the distribution 
of the current density j would depend not only on the 
radial coordinate p but also vary along the poloidal 6 and 
the toroidal ip angles due to the presence of the ( m/n = 
1/1) magnetic island. This agrees with the analysis of 
numerous disruptions in the JET tokamak J25| . One can 
assume that the radial profiles of the RE current density 
averaged along poloidal and toroidal angles are almost 
uniform. This gives the value of the safety factor at the 
beam axis g(0) is less than unity. This assumption is 
supported by a number of experimental measurements 
of the current profile after the sawtooth crashes in the 
TEXTOR, the TFTR, and JET tokamaks . 

The toroidal electric field accelerates electrons to 
higher energies. With increasing electrons energy their 
orbits drifts outwardly [H, [33| and eventually hit the 
wall. It is illustrated in Fig. Q] (a). This effect may be 
one of mechanisms of slow RE current decay. Calcula¬ 
tions show that the outward drift velocity Vdr is of the 
order of a few m/s for typical discharges in TEXTOR. 
The RE current decay rate dl p /dt due to outward drift 
RE orbits can be roughly estimated as follows. This loss 
mechanism is mainly caused by the shrinkage of the beam 
radius a. The rate of such a shrinkage da/dt is of the or¬ 
der of the average outward velocity Vdr■ Since I p oc a 2 , 



FIG. 4: (a) Evolution of a RE orbit in the (R, Z)-plane for 
the TEXTOR discharge # 117527. (b) Poincare section of 
RE orbit of energy 11.7 MeV in the ( R , Z)-p\ane. The plasma 
current I p = 50 kA. 


we have dl p /dt oc (2 I p /a)da/dt = (2 I p /a)vd r - For the 
typical values of I p ~ 0.2 MA, a ~ 0.2 m, and Vdr ~ 1 
m/s one has dlp/dt ~ 4 MA/s. This estimation is in the 
order of the experimentally measured average decay rate 
of the runaway current plotted in Fig. [2] (b). 

The effect of magnetic perturbation on RE beams de¬ 
pends on their safety factor profile q. The latter varies 
in the interval [g(0) < 1 ,q{a)\ with its edge value q(a) 
less than 3/2 [or 4/3, 5/3]. Such a RE beam is relatively 
stable to the effect of magnetic perturbations. The single 
m/n = 1/1 mode does not create the stochastic layer at 
the beam edge for REs with energies up to several MeVs 
since their drift surfaces are close to magnetic surfaces. 
With increasing the energy of electrons the drift surfaces 
strongly deviates from magnetic ones and thus creates 
the perturbation harmonics with higher mode numbers 
m > 1. The interactions of several resonance modes of 
perturbations may form the stochastic zone at the beam 
edge which leads to fast RE losses as illustrated in Fig. 0] 
(b). This process may explain the sudden RE current 
drop accompanied by magnetic activity and RE bursts 
observed in experiments (see, e.g., [j,jj|). 

Summary. Based on the analysis of numerous experi¬ 
mental data obtained in the TEXTOR tokamak we have 
proposed the mechanism of RE beam formation during 
the plasma disruption. The plasma disruption starts due 
to a large-scale magnetic stochasticity caused by nonlin- 
early excited of MHD modes with low (to, n) numbers 
{m/n = 1/1, 2/1, 3/2, 5/2,...). The RE beam is formed 
in the central plasma region confined by the intact mag¬ 
netic surface. Its location depends on the safety factor 
profile q{p) and the spectrum of MHD modes. In the 
cases of plasmas with the monotonic profile of q{p) and 
at sufficiently small amplitude of the m/n = 1/1 mode 
the most stable RE beams are formed by the intact mag¬ 
netic surface located between the magnetic surface q = 1 
and the closest low-order rational surface q = m/n > 1 
{q = 5/4, q = 4/3, or q = 3/2). 

This mechanism reproduces well the essential features 
of the measurements. Particularly, the TQ and the CQ 
are determined by the strong electron diffusion and am¬ 
bipolar transport of particles in a stochastic magnetic 
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field, respectively. The slow decay of the RE current 
is due to the outward drift of RE orbits induced by a 
toroidal electric field, and the spiky quick decay of REs is 
due to resonant interaction of high-energy REs with the 
m/n = 1/1 MHD mode. The effect of external resonant 
magnetic perturbations on low-energy electrons (up to 5- 
10 MeV) is weak and does not cause their loss. This is in 
agreement with the recent experiments in the TEXTOR 
tokamak §1] . The detailed description of the mechanism 


of RE formation and the evolution of RE current based 
on the analyses of experimental observations will be given 
in a separate publication [35ij. 
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